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Abstract 

The angular distributions of the unpolarised differential cross section and tensor 
analysing power A^x of the dd ^ ar] reaction have been measured at an excess 
energy of 16.6 MeV. The ambiguities in the partial-wave description of these data 
are made explicit by using the invariant amplitude decomposition. This allows the 
magnitude of the s-wave amplitude to be extracted and compared with results 
published at lower energies. In this way, firmer bounds could be obtained on the 



Preprint submitted to Nuclear Physics A 



20 November 2008 



scattering length of the r]a system. The results do not, however, unambiguously 
prove the existence of a quasi-bound t] a state. 
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1 Introduction 

The possibility of the i] meson forming a quasi-bound state in a nucleus was 
first raised by Haider and Liu [1]. Such a state could arise as a consequence of 
the strongly attractive ?7-nucleon interaction that is driven by the A^*(1535)S'ii 
resonance. Since the rj has isospin / = 0, this leads to attraction for both the 
protons and neutrons in the nucleus. A value of the (complex) s-wave rjN 
scattering length a^Tv = dr +icii ~ (0.28 + 0.19i) fm had been found [2] on the 
basis of the existing phase shifts. Using this value, it was shown that the r] 
meson could form a quasi-bound state with nuclei of mass number A> 10 [1]. 
Other groups found similar results when starting from this relatively small 
value of ttrjN [3,4]. 

The likelihood of ?7-nucleus quasi-bound states existing would increase signifi- 
cantly if the ?7-nucleon scattering length were much larger than that assumed 
by Haider and Liu in 1986. In the subsequent years, widely differing estimates 
have been given for o^at, with some real parts being as large as 1 fm; see Ref. [5] 
for a summary. Theoretical studies by Ueda [6] gave some hints that the t] me- 
son might form a quasi-bound state with ^'^He and even with deuteron. The 
states could be narrow and lightly bound in few-nucleon systems and therefore 
might be observed through their effects above threshold. 

According to the Watson-Migdal theory [7,8], when there is a weak transition 
to a system where there is a strong final state interaction {FSI), one can 
factorise the s-wave reaction amplitude, fs, near threshold in the form 



where is the rj cm. momentum. The unperturbed production amplitude 
Jb is assumed to be slowly varying and is often taken to be constant in the 
near-threshold region. 
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Unitarity demands that the imaginary part of the scattering length be positive, 
i.e., tti > 0. In addition, to have binding, there must be a pole in the negative 
energy half-plane, which requires that [5] 



ai\/ \ar\ < 1 . 



(1.2) 



Finally, in order that the pole lie on the bound- rather than the virtual-state 
plane, one needs also < 0. 

The unexpectedly large near-threshold production amplitude in the case of 
the pd — >■ r^^He reaction, as well as its rapid decrease with rising energy [9, 
10], were interpreted as evidence for a strong s-wave FSI, which might be 
associated with the formation of a 77 ^He quasi-bound system [11]. Analogous 
experimental data on the dd — > 77 ^He reaction [12, 13] show an amplitude that 
varied more slowly than in the case of 77 ^He. A combined analysis to the two 
data sets within the framework of a simple optical potential yielded scattering 
lengths [13] 



Taken at face value, these suggest that there are poles in both amplitudes 
close to threshold but that, whereas the 77 ^He system might be quasi-bound, 
the scattering length for 7/^He does not satisfy the condition of Eq. fll.2p . 
Nevertheless, more recent very refined measurements of the pd — rj^He dif- 
ferential cross section show even larger 7/^He scattering lengths [14,15]. These 
prove that there is a pole in the complex energy plane within about 1 MeV 
from threshold and this is confirmed by the energy dependence of the angular 
distribution [16]. 

In order to isolate the s-wave amplitude for dd — > 7/^He, it is necessary to 
measure differential and not merely total cross sections [12,13]. This has been 
done by the ANKE collaboration at an excess energy of Q = 7.7 MeV [17]. 
However, as shown there, even differential cross sections are not completely 
sufficient for this purpose and the analysis of these and the earlier data is 
model dependent. Although there is only one s-wave amplitude, observed de- 
viations from isotropy might be a purely p-wave effect or could result from an 
interference of the s- with the rf-wave [17]. The only way to decide between 
these two solutions is through the determination of a deuteron analysing power 
in conjunction with the differential cross section. 

In this paper we present results of a measurement of da/dQ and the tensor 
analysing power A^x a.t Q = 16.6 MeV, which allows us to resolve this ambi- 
guity. The experiment was carried out using the Big Karl spectrograph [18,19] 
and the basic setup is described in section It was carried out in two stages. 



a(r/^He) 
a(r/^He) 



(-2.3 + 3.2i) fm, 
(-2.2 + l.lz) fm. 



(1.3) 



3 



first using an unpolarised beam, as discussed in section [2l2l In order to extract 
an analysing power, the beam polarisation had to be established and the way 
of doing this is explained in section 12.31 and then used in the measurement of 
A^x in section 12.41 

The results for the unpolarised differential cross section are presented in sec- 
tion l3.1l and analysed in terms of the partial wave and invariant amplitudes [17] 
in section [3^ It is shown that, if we neglect g-waves in the rja system, the 
magnitude of the s-wave amplitude could be determined from measurements 
of differential cross section and A^^. The data clearly favour s-d interference 
as the origin of the angular dependence. The results are compared in sec- 
tion 13.31 with those obtained at lower energies and together these go a long 
way in clarifying the ambiguities regarding the s-wave scattering length. Our 
conclusions are drawn in section |H 



2 Experiment 

2.1 Experimental Setup 

The GEM collaboration measured the dd ^ arj reaction at the COSY-Jiilich 
accelerator using vector and tensor polarised deuteron beams. The beam mo- 
mentum of 2385.5 MeV/c corresponds to an excess energy oiQ = 16.6 MeV for 
this reaction when an //-meson mass of = 547.7 MeV/c^ is used [20]. The 
Big Karl magnetic spectrograph [18, 19] employed for this study is equipped 
with two sets of multi-wire drift chambers (MWDC) for position measure- 
ment and thus track reconstruction. Two layers of scintillating hodoscopes, 
4.5 m apart, led to a more accurate time-of-fiight measurement than previ- 
ously achieved with Big Karl. They also provided the energy-loss information 
necessary for particle identification. 

A sketch of the target area is shown in Fig. [TJ The target was a cell of 6 mm 
diameter and 4 mm length, which contained liquid deuterium or hydrogen. The 
target windows were made from Mylar with a thickness of less than 1 fim [21]. 
Crossed scintillator paddles, working in coincidence, were used as luminosity 
monitors. They were calibrated by using scattered particles from a beam with 
reduced intensity that were directly counted in an exit window in the side 
yoke of the dipole Dl on the way to the beam dump. The detectors were 
therefore named the left and right beam-dump monitors. Whereas the beam 
telescope saturates for intensities above about 10 ^s~^, the luminosity monitors 
were still far below saturation at beam intensities of the order of 10 ^s~^. 
The luminosities obtained from the two monitors always agreed within a few 
percent. 
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M detector 



Fig. 1. Schematic layout of the experimental set-up close to the target area. Left: 
Positioning of the luminosity counters. Right: View of the M detector counter used 
for monitoring the beam polarisation. 

The system yielded very good particle separation, as is seen from the two- 
dimensional plot of energy loss versus time of flight shown in Fig. [2l 
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Fig. 2. Energy losses in the scintillator layers versus the time of flight between 
them for the d{d,X) reaction at 2.39 GeV/c. The islands of "^He and ^He are well 
separated from the recoil particles. 



2.2 Unpolarised Angular Distribution 



In the first stage of the experiment, we used an unpolarised deuteron beam 
incident on the target cell filled with liquid deuterium. The acceptance of the 
Big Karl spectrograph for a-particles is limited, not only in momentum, but 
also geometrically, mainly due to the side yokes in the first dipole magnet. The 
consequences of this for the acceptance are shown in Fig. [3] in terms of the 
polar {9) and azimuthal (0) angles of the a-particles. The regions with high 
acceptance that were retained in the analysis, tt/2 — 0.5 < |0| < tt/2 + 0.5, 
are indicated by the solid lines. The simulated acceptance with these cuts is 
a smooth function ranging from 0.3 at cos 6' = 0.4 to 0.9 at cos 6' = 1. 
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Fig. 3. Measured a-particle events after applying cuts as function of the polar angle 
and the cosine of the azimuthal angle in the center-of-mass frame. The regions 
around (/) ~ and (p « ztyr have almost zero acceptance due to the side yoke in the 
first magnetic dipole. 

After applying cuts, the missing mass of the detected a-particle was evaluated 
and the event placed in a cos 6 bin. The resulting spectrum, summed over all 
angles is, shown in Fig. |H A clear t] peak is superimposed on a background 
coming from multipion production. The latter was simulated by calculating 
the spectral form of two, three and four pion production in a phase-space 
model, though the shapes from two and three pion production were almost 
identical. The individual strengths were fitted to the experimental spectrum 
together with a Gaussian for the peak and the results shown in Fig. HI 

jo= 

zr ++++++++++ experimental data for d+d — > "He+X 
background simulation 
background+gaussian simulation 
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Fig. 4. Missing mass distribution for the dd — > aX reaction in the range 
0.4 < cose* < 1. 

An analogous analysis was applied to six bins, each of width A(cos^) = 0.1. 
The peak areas were corrected for acceptance and converted into differen- 
tial cross sections using the known luminosity. The resulting differential cross 
sections are shown in Fig. [51 

When the incident deuterons are unpolarised, the entrance channel is sym- 
metric and a representation of the differential cross section in terms of cos^ 6 
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Fig. 5. Unpolarised dd ^ arj differential cross section. Since tlie entrance cliannel is 
symmetric, the data are shown in terms of cos^ 6. Upper panel: Data from Ref. [17] 
taken at Q = 7.7 MeV. The curve is a linear fit in cos^ 6. Lower panel: Same as the 
above but for the present data at Q = 16.6 MeV. These data are compared with a 
quadratic fit in cos^ 6. 

is appropriate. Already at an excess energy of 7.7 MeV the ANKE collabo- 
ration found deviations from isotropy [17], with evidence of a forward dip. In 
contrast, the present data at Q = 16.6 MeV show a strong peaking in the 
forward direction. Both sets of results are shown in Fig. [5l We defer further 
discussion of these distributions until section [3l 



2.3 Beam Polarisation 

The ion source used to generate the polarised deuteron beam is described 
in Ref. [22]. For this experiment, in addition to the unpolarised beam, two 
different combinations of vector pz and tensor pzz polarisation were used ^ I: 



It is customary to call the quantisation axis "z" in the source frame. 
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Pz = -1/3 and p^^ = +1, 
p^ = -1/3 and p^^ = -1. 



(2.1) 



The vector polarisation was measured with a low energy polarimeter placed 
in the injection beam line, where the deuteron energy is about 76 MeV. Using 
a carbon target, deuteron elastic scattering was measured with scintillating 
detectors. The results for the two polarisation combinations were consistent 
with the expectations, with = —0.32 ± 0.02 and Pz = —0.33 ± 0.02. It has 
been shown that there is little or no loss of polarisation in the subsequent 
acceleration of deuterons at COSY up to the momentum of interest here [23]. 

The tensor polarisation of the beam was measured by scattering the accel- 
erated deuterons from the liquid hydrogen target. The elastically scattered 
deuterons were then measured with Big Karl. The cross section for polarised 
dp elastic scattering in the backward direction can be expressed as: 




180°) 



pol 




180°) 



unpol 



2Pzz^yy{ 



180°] 



.(2.2) 



At strictly 180°, there is only one analysing power, Ayy = A^^, and this was 
measured at Saclay [24] over a broad range of deuteron momenta that covered 
the 2.39 GeV/c used here. The points were remeasured with smaller uncer- 
tainties by Punjabi et al. [25], and we made use of their results. 
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Fig. 6. Missing mass spectrum for the dp — > pd reaction for the unpolarised beam 
(solid line) and two different tensor polarisations (dashed and dotted lines). 

The outgoing deuterons were identified using the specific energy losses and 
time of flight. In order to eliminate background, the deuteron missing-mass 
spectra were analysed. An example of such a spectrum is shown in Fig. [6] for 
the dp — s> dX process for the unpolarised beam and for both states of deuteron 
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polarisation. In all cases there is a strong proton peak sitting on a smoothly 
varying background from which it was straightforward to extract the num- 
bers of elastically scattered deuterons. The tensor polarisations summarised 
in Table [T] were then determined on the basis of Eq. (12.21) . 

Table 1 

Results of the measurements of vector and tensor polarisations of the deuteron 
beam. The first errors are statistical and the second systematic, arising principally 
from the calibration polarisations [24, 25] . 



Pz 


Pzz 


nominal 


measured 


nominal 


measured 


-1/3 


-0.33 ±0.02 


-1 


-0.87 ±0.11 ±0.01 


-1/3 


-0.32 ±0.02 


±1 


+0.91 ±0.14 ±0.01 



2.4 Polarised Cross Section Measurement 



In the subsequent experiment, we employed the vector and tensor polarised 
deuteron beam incident on the deuteron target to measure the analysing power 
in the dd ^ at] reaction. The following sequence of polarisation modes was 
used 

V+iV-iV+iV^iV+iV-iV+iV-iV^V+iV-i • • • , (2.3) 

where "Pj, 2=+!, —1, is a mode with the nominal values of vector and 
tensor p^z polarisation of the beam given in Table [H It should be noted that, 
in contrast to the initial experiment with the unpolarised beam, the nominally 
unpolarised mode was prepared by employing three transitions of the source 
and so some small residual polarisation could not be excluded. 

In the reaction frame, z is taken along the beam direction, y represents the 
upward normal to the COSY accelerator, and x = y x z. The polarised dd — 
a rj cross section at a production angle 6 can be written in Cartesian basis as 



da ,„ , A I f do 



dVl J pol / V dVl J unpol 



%{d)) =[l±^A,(%.COS0 



2' 

H 



4 ^ ^ ' 4 ^ ' 



(2.4) 



G 



The polarisation was inspected for each run by fitting the quantities F, G, 
H and J to the angular distributions on the M detector. An example of the 
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response of the M detector to different beam polarisations is shown in Fig. [71 
together with the fitted curves. 
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Fig. 7. Ratio of the counting rate in each M detector segment, corresponding to a 
certain azimuthal angle (p, for two different polarisations of the deuteron beam. The 
dashed line represents the fit result for {pz,Pzz) = (— while the solid line 
corresponds to {pz,Pzz) = (~|) 
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Fig. 8. Comparison of the ratio of the measured parameters H and J to the reference 
values i^i-ef and Jrcf for both polarisation states as functions of the run number. The 
H parameter corresponds to the vector polarisation pz, while J is connected with 
the tensor beam polarisation pzz- Left: Polarisation state with pz = — |, Pzz = +1- 
Right: Polarisation state with pz = — ^, Pzz = — 1- 

The fitted parameters H and J are shown in Fig. [H] relative to the nominal 
values i^ref and Jref for individual run numbers. The deviations from these 
values are less than ±15%. Because of the smallness of these fluctuations, we 
retain the nominal values in the analysis. 

The analysing powers could, in principle, be extracted by fitting Eq. fl2.4l) di- 
rectly to the data. However, the present detector does not have full acceptance 
(see Fig. [3]). Furthermore, as is seen from Eq. (12.31) . the unpolarised beam in 
this part of the experiment had only 10% of the statistics and even here there 
may be some residual polarisation. We therefore developed a method to ex- 
tract the analysing powers without making use of the unpolarised beam. 



10 



As indicated in Fig. [3l there are two ranges of azimuthal angle where there is 
almost complete acceptance, i.e., 



i(7r- 1) < 101 < i(7r + l)rad. (2.5) 

For these intervals the mean value of cos cf) vanishes and that of cos 20 is given 
by 

<cos20>=/ cos20rf0 = -0.84. (2.6) 

J(7r-l)/2 

It therefore becomes clear that Ay could not be measured with the present 
layout. 

Since < cos 20 > is close to —1, Eq. (12. 4p shows that our experiment is primar- 
ily sensitive to the values of A^^, with only a small contamination from Ayy. 
We therefore neglect the contribution of A^x — Ayy to the counting rates in 
these intervals. In this approximation, the cross section integrated over these 
intervals of azimuthal angle becomes simply 

where the unpolarised cross section is integrated over the same range. 
Carrying out this procedure for the two polarisation states, we find that 

^ ^ - Ii ^ 0.23 (P+ .2 8) 

/+ + /- l + 0.23A,,(p+ +p-) 

Using the measured beam polarisations reported in Table [H we can solve to 
find the value of A^x'- 

Axx = 2.44 A/(l - 0.02 A). (2.9) 



Since the unpolarised cross section drops out of Eq. (12. 9p . it is sufficient to 
count the numbers of registered a t] events per incident beam n. The total 
numbers of beam particles for the truly unpolarised beam and the two polar- 
isation states are given in Table [2l Uncertainties in the target thickness and 
in the acceptance etc. cancel. 

As is demonstrated later, A^x as well as the unpolarised cross section are 
both even functions of cos 6 and so they are shown in Fig. [9] as functions of 
cos^ 6. The large error bars for Axx arise from the fluctuations in the significant 
background under the rj peaks. It should be noted that the analysing power 
Axx in the range just above the rj peak is found to be consistent with zero. 
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Table 2 

Integrated beam intensities for the unpolarised and polarised deuteron beams. 



Nominal 

Pz Pzz 


Integrated beam 
intensity n 




4 +1 
-I -1 


(12.70 ±0.48) X 10^3 
(6.98 ±0.62) X 10^3 
(6.19 ±0.51) X 10^3 



< 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
COS^O) 



Fig. 9. Upper panel: Differential cross section for the dd ^ arj reaction. Lower panel: 
Analysing power Axx- The solid curves represent a fit with four partial waves; the 
dotted curves with invariant amplitudes. 

3 Results and Discussion 

3.1 Unpolarised Cross sections 

In order to deduce the total dd —>■ ar] cross section, the differential data were 
fitted in terms of Legendre polynomials: 

-7^= Y^aMcose). (3.1) 

The smallest value for 2imax required to describe our data was found to be four. 
The resulting parameters are Oq = 1.27 ±0.03 nb/sr, 02 = —0.29 ±0.06 nb/sr, 
and a4 = 1.65 ± 0.07 nb/sr, and the corresponding fit is shown in Fig. [51 
This indicates that there must be at least d-wave contributions. This is to be 
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contrasted to the lower energy ANKE results [17], where 2irnax = 2 suffices 
(see Fig. [5]). In their case, oq = 1.30 ± 0.18 nb/sr and 02 = —0.79 ± 0.19 nb/sr 
which, as in our results, is negative. 

The total cross section deduced from the fit is 



a = 16.0 ± 0.4 nb, 



(3.2) 



where uncertainties in the target thickness, incident flux, and acceptance in- 
troduce an additional systematic error of ±1.6 nb. In Fig. [10] this result is 
compared with existing data, which are all taken at smaller excess energies. 
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Fig. 10. Excitation function for the total cross section for the dd ^ arj reaction. 
The data are from Frascaria et al. [12], Willis et al. [13], and Wronska et al. [17]. 
Only statistical errors are shown. 



Willis et al. [13] only extracted the cross section for helicity m = ±1, which 
makes their result for the unpolarised total cross section model dependent. 
The values shown in Fig. [10] correspond to the assumption that the p-wave 
amplitude leading to the m = cross section is negligible. This is consistent 
with our data and we will return to the point in section 13. 3[ after discussing 
the amplitude decomposition. 



3.2 Amplitudes and Observables 



In order to extract the s-wave amplitude for dd ^ ar], we attempted to fit 
partial wave amplitudes 

ai = \a,\e"^\ (3.3) 
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to the present data[l]- Since the angular distributions require terms up to 
cos"^ 9, we consider s, p and (two) d waves. After applying angular momentum 
and parity conservation, and taking into account the identical nature of the 
incident deuterons, it is seen that there are only four transitions, which are 
noted in Table [31 

Table 3 

Quantum numbers corresponding to the four lowest partial waves for the dd — > 
ar/ reaction. Here Si is the total spin of the initial deuterons with orbital angular 
momentum ^i. Bose symmetry requires that Sj + £j be even. The total angular 
momentum J is equal to the a rj orbital angular momentum i f . 



amplitude name 


Si 






wave name 


ao 


1 


1 





s 


ai 


2 


2 


1 


P 


a2 


1 


1 


2 


d 


03 


1 


3 


2 


d 



Simonius has shown how to relate such amplitudes to different possible ob- 
servables [26]. However, when these are fitted directly to our data it is found 
that only the magnitude of the s-wave amplitude |ao| is stable. In contrast, 
the magnitudes of the two d-wave amplitudes |a2| and jasl are completely cor- 
related. To see the origin of these effects, it is simplest to consider the relation 
of the partial waves to the invariant amplitudes used by Wrohska et al. [17], 
which we briefly summarise here. 

Due to the identical nature of the incident deuterons, only three independent 
scalar amplitudes are necessary to describe the spin dependence of the reac- 
tion. If we let the incident deuteron cms momentum be pd and that of the rj 
be p^, then one choice for the structure of the transition matrix M. is 



M = A{ei X €2) ■ pd + B{ei x €2) ■ [pd x (p^ x pd)] {pr, ■ Pd) 

+C [(el • Pd) 62 ■ {Pr, X Pd) + (62 ■ pd) Ci ■ {prj X p^)] , (3.4) 

where the q are the polarisation vectors of the two deuterons. 

If we retain up to d waves in the final system, then B and C have no angular 
dependence while that of the amplitude A can be written as 

A{e) = + A2P2{co^e). (3.5) 



^ These amplitudes should not be confused with the scattering length, also denoted 
by a. 
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The observables in a spherical basis were related to these amplitudes in Ref. [17] 
but they can be easily converted to the Cartesian observables used here [27]. 
Our experiment only yielded data on the unpolarised cross section and tensor 
analysing power A^^^ and these are given in terms of the invariant amplitudes 
by 



7^ = - d^ol' + 2/?e(AoA;)P2(cos 9) + 1^2^ (Palcos^))' 



dfi Pd 



[I + ^ = 2^(\B\' sin2 6 cos^ + \C\^ sin' 6 



(3.6) 
(3.7) 



where the results have been expressed in terms of convenient linear combi- 
nations. From these it is seen that both the cross section and A^x are even 
functions of cos 6*, a result that has already been used in our analysis and 
presentation. The data could therefore in principle fix the magnitudes of the 
amplitudes Aq, A2, B, and C, and the interference between Aq and A2, while 
being completely insensitive to all the other phases. Furthermore, the linear 
combinations of Eqs. (13. 6p and (13. 7p show that the fitting of |-B| and |C| is de- 
coupled from that of A. The parameters resulting from fitting the data in this 
basis are given in Table HJ with the fit curves being shown in Fig. O From this 
it is seen that the A amplitudes are dominant, with C being consistent with 
zero within error bars. If |-B| also vanished, it would follow from Eq. (13. 7p that 
Axx = ~| for all angles. On the other hand, even the small contribution from 
the B term changes the angular dependence of Axx, as is evident in Fig. [9l 



Table 4 

Fit results of invariant amplitudes of Eqs. (j3.6p and (|3.7p to the present data. Since 
\C\ was found to be zero within error bars, it was put exactly to zero. 



fit parameter 


value 


\Ao\' 


6.6 ± 1.7 


2Re {A*qA2) 


-25.0 ±9.5 


\A2\' 


48.4 ± 14.5 




9.3 ±5.1 


\c\' 






The Wrohska et al. amplitudes can be easily expressed in terms of those of 
the partial waves of Table [3l Explicitly 
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From the above discussion, it is seen that if £^ = 4 waves are neglected, then 
the data might determine the absolute magnitudes of the = and = 1 
amplitude, viz. \ao\ and |ai|. On the other hand, only two linear combinations 
of the d-wave amplitudes |a2 + 203/ a/14 P and \a2 — 3a3/A/l4p could fixed 
by the data. Thus the differential cross section and tensor analysing power 
Axx are together only sensitive to 7|a2p + 3|a3p rather than the individual 
magnitudes. This accounts for the complete correlation found in the direct 
fitting of the partial wave amplitudes to the data. 



3.3 The s-Wave Amplitude 



We have seen that, provided waves are neglected, the present data allow us 
to extract the magnitude of the s-wave amplitude |ao|. From this we obtain a 
spin-averaged square of the s-wave amplitude, l/^p through 

dcTs _P'n\f \2 _ 2^7? I 4 |2 _ 1 1 I 12 (nQ\ 

dil Pd 3Pd 2747r 

Using the value given in Table HI we find that l/^p = 4.4 ±1.1 nb/sr. 

Close to threshold, where s-waves dominate, l/^p can be extracted from the 
total cross section a via 

l/.P = -f, (3.10) 
p^ Att 

since the contribution from s-d interference drops out. 

Equation fl3.10p has been used by Willis et al. [13] to extract but there 
are two important considerations here. Our data show that the B and C 
amplitudes play only a minor role in the observables and this is likely to be 
even more true at lower energies. Hence to a good approximation the m = 
cross section can be neglected and the total cross section should be two thirds 
of their m = ±1 value. This result has already been used when plotting the 
data in Fig. [5l 
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Furthermore, although the SPESIII detector had complete acceptance for the 
a particles from rj production, the full polar angle could not be reconstructed 
and the authors assumed isotropic distributions [13]. While this is a good 
approximation close to threshold, it is not valid at their two highest energies, 
as is evident when comparing with the data of Ref. [17]. There will therefore be 
(i-wave contributions which must be subtracted in Eq. (13.91) before extracting 
l/sp. To do this, we assume that the d-wave amplitudes A2 and B vary with 
a threshold factor of p^, where the normalisation is fixed by our results given 
in Table H 



This procedure results l/^P = 13.8 ± 1.2 nb/sr at = 73 MeV/c and 10.6 ± 
1.3 nb/sr at pr, = 91 MeV/c for the Willis data [13] and 14.3 ± 2.4 nb/sr for 
the Wrohska measurement at 86MeV/c [17]. All the values of l/^P are then 
shown in Fig. [TTl 



Willis et al. [13] did a combined optical-model fit to all the near-threshold 
pd — ^ ?7'^He and dd — 77 ^He data as discussed in the introduction, assuming 
that only s-wave production occurred, and found a value of the 77 a scattering 
length of a = (—2.2 + il.l) fm. Given the extended dd — > 77 ^He data set now 
available, we have attempted a direct fit in the scattering-length approxima- 
tion 



|f|2 = (311) 



The best fit shown in Fig. [TT]is obtained for \ar\ = 3.1 ±0.5 fm, = 0±0.5 fm 
and = 34 ± 1 nb/sr, corresponding to a quasi-bound or virtual state with 
\Qo\ f» 4 MeV. These above-threshold data are, of course, insensitive to the 
sign of ttr so that they could never tell whether the system is quasi-bound or 
virtual. The argument given in Ref. [13] is that, since 77 ^He is more likely to be 
bound than r/^He, the fact that IQol is smaller for 7/^He [10,14,15], suggests 
that T^'^He is indeed quasi-bound. 



To try to explore the systematic uncertainties, one can include also an effective 
range term r, as has been done for the near-threshold dp — > 77 '^He data [14] 
and fit 

fs = , ■ \i 2- (3-12) 

However, due to the strong coupling between the fit parameters, the error bars 
become exceedingly large, = 6.2 ± 1.9 fm and Oj = 0.0.001 ± 6.5 fm. This fit 
is also shown in Fig. [11] The large errors here do not allow a decisive answer 
on the position of the r] a pole. 
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Fig. 11. The world data for spin- averaged square of the magnitude the s-wave am- 
phtude for dd — > r/^He as a function of the r] cm momentum. The dashed curve 
shows the best fit Eq. ()3.1ip with parameters given in the text. The sohd curve is 
a fit with Eq. (|3T2D . 



4 Summary 



We have measured the differential cross section and tensor analysing power 
Axx of the dd ^ arj reaction at an excess energy of 16.6 MeV. The recoil- 
ing «-particles were measured in a magnetic spectrograph and the 77 mesons 
identified through a missing-mass technique. The biggest uncertainty in the 
method arises from the large background coming from multipion production. 
Despite this, an angular distribution and a total cross section could be given. 

The angular distribution of the analysing power A^x was measured with a 
tensor polarised beam. Since this observable required measurements with two 
different spin modes, the uncertainty due to the background subtraction is even 
larger. Nevertheless, the values thus obtained showed that p-wave amphtude 
was very small so that deviations from isotropy must come primarily from s-d 
interference [13,17]. 

Combining these results with the unpolarised cross section data allowed the 
spin-averaged square of the s-wave amplitude to be extracted with rea- 
sonable error bars. Assuming that the d-wave production cross section varied 
like Q^, values of |/<jp were also deduced from earlier measurements such that 
there are now 13 values for Q < 16.6 MeV. Fitting the momentum depen- 
dence in the scattering length approximation [7,8] gives \aj.\ = 3.1 ± 0.5 fm 
and Oj = 0.0 ± 0.5 fm. This suggests that there is a singularity in the complex 
plane close to threshold but its position is far from clear and the uncertainty 
grows if a fit is attempted with an effective range as well as a scattering length. 
Further data are clearly needed! 
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